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ABSTRACT 

The  subject  of  this  paper  is  the  buckling  behavior  of  a  rectangular  plate,  with  parallel  thin- 
walled  stiffeners  attached  to  one  side,  subjected  to  a  combination  of  axial  compression,  lateral 
pressure,  and  bending  moment.  The  plate  is  modeled  by  the  Von  Kdrmdn  plate  equations  and 
the  stiffeners  by  a  nonlinear  beam  theory  recently  derived.  An  analytical  solution  is  obtained 
for  the  buckling  load  corresponding  to  a  torsional  tripping  mode  of  the  stiffeners.  The  effects 
of  various  boundary  conditions,  imperfections,  and  residual  stress  are  included. 

1.     INTRODUCTION 

Stiffened  plates  are  a  basic  structural  component  of  ships  and  submarines.  These  struc- 
tures are  designed  with  generous  safety  margins  against  overall  collapse  triggered  by  buckling. 
The  object  of  analytical  work  is  to  determine  design  criteria  to  inhibit  buckling  at  any  stress 
less  than  yield.  Recently  [see  Danielson  et  al  (1993)],  we  have  developed  an  analytical  formula 
for  the  buckling  load  of  a  stiffened  plate  subjected  to  a  combination  of  axial  compression 
and  lateral  pressure.  The  object  of  the  present  paper  is  to  improve  and  extend  our  previous 
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analysis.  A  review  of  the  literature,  given  in  our  earlier  work,  will  not  be  repeated  here  but 
details  of  the  analysis,  which  supersedes  our  earlier  work,  will  be  recorded  here. 

We  first  consider  a  plate  which  is  initially  rectangular  in  shape  and  has  several  parallel 
I-stiffeners  spaced  a  distance  6  apart.  The  structure  is  subjected  to  a  combination  of  uniform 
axial  compressive  stress  a  (force  per  unit  area  of  a  side),  uniform  lateral  pressure  p  (force  per 
unit  lateral  area  of  the  plate),  and  uniform  bending  moment  M  (moment  per  unit  length  of 
an  edge).  We  suppose  that  at  low  values  of  <r,  p,  and  M  the  plate  and  stiffeners  simply  bend 
and  compress  symmetrically.  Our  object  is  to  find  the  critical  load  at  which  the  stiffened 
plate  may  buckle  into  an  alternate  mode  (see  Fig.  1). 

Our  present  analysis  is  based  on  the  following  simplifying  assumptions: 

(i)  Each  plate-stiffener  unit  of  width  b  undergoes  an  identical  deformation. 

(ii)  The  plate  obeys  the  nonlinear  Von  Karman  plate  equations  [see  Timoshenko  and  Gere 
(1961)].  The  stiffeners  obey  the  nonlinear  beam  equations  derived  by  Danielson  and 
Hodges  (1988). 

(iii)  The  plate  and  stiffener  material  is  elastic,  linear,  and  isotropic. 

(iv)  Every  particle  on  the  bottom  surface  of  a  beam  undergoes  the  same  displacement  as 
the  corresponding  particle  on  the  top  surface  of  the  plate,  and  every  line  of  particles 
in  the  beams  normal  to  the  plate  surface  remains  normal  to  the  deformed  plate  at  its 
surface.  In  other  words,  the  bases  of  the  stiffeners  are  clamped  to  the  plate. 

(v)  The  prebuckling  displacements  are  less  than  the  maximum  thickness  of  the  structure 
and  independent  of  the  transverse  coordinate. 

(vi)  The  incremental  buckling  extensional  strains  at  the  midsurface  of  the  plate  are  negli- 
gible. 


(vii)  The  incremental  buckling  displacements  may  be  approximated  by  the  fundamental 
harmonic  in  their  Fourier  expansions. 

(viii)  The  plate  and  beams  are  so  thin  that  their  thicknesses  are  negligible  compared  to 
their  width,  height,  length,  and  the  wavelength  of  deformation.  A  stiffener  is  so  slender 
that  its  width  and  height  are  negligible  compared  to  its  length  and  the  wavelength  of 
deformation. 

2.  GENERAL  POTENTIAL  ENERGY  FUNCTIONAL 

It  follows  from  assumption  (i)  that  we  need  only  analyze  a  single  plate  unit  containing 
a  single  stiffener.  From  assumptions  (ii)-(iii),  the  potential  energy  of  the  plate  plus  beam  is 
given  by: 

■*    «-    Et     -2  -2 


+  atui  —  pw  +  Mwu  >  dxi  dx2 


\w2  w2 

+D  Mi  +  uwuw22  +  -2-  +  (1  -  v)w22 
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(1) 
Here  (xj ,  x2,  x3)  are  Cartesian  coordinates  measured  from  the  midpoint  of  a  side  of  the  plate. 

The  plate  unit  has  length  a,  width  6,  and  thickness  t,  while  the  beam  has  cross-sectional 

area  A  and  centroidal  height  c.  The  elastic  constants  are  defined  by 


12(1  -z/2)1  2(1  +  //) 

where  E  is  the  Young  modulus  and  u  is  the  Poisson  ratio.  The  displacements  of  the  plate 
midsurface  in  the  X\,x2,  and  23  directions  are  denoted  by  u(xi,  x2),  v{x\,x2),  and  iu(ii,  x2), 
respectively.  Subscripts  on  u,v,  or  w  denote  partial  differentiation  with  respect  to  the 
coordinates  X]  or  x2;  e.g.,  wi2  =  dx  qt  ■  The  extensional  strains  at  the  midsurface  of  the 
plate  are  given  by 


1      2  1/  x  1      2 

en  =  ui  +  -iflj,  ew  =  -(W2  +  fi  +  wiw2),  e22  =  u2  +  ^2 

The  strains  in  the  beam  are  denoted  by  7n(xi,  x2,  *3)»7w(*ii  x2,  x3),  and  713(1!,  x2,  x3) 
and  are  related  to  the  displacements  by  equations  (9,  10)  of  Danielson  and  Hodges  (1988), 
which  upon  invoking  assumption  (iv)  are  transformed  into: 

7n     =    Eu  +  (Eu  +  e"i2)<£3  -  £13^2  +  2^2  +  9^3 
7l2     =     £12  -  yEU(f>3,  713  =  £13  +  -Eufa 

Eu     =     en  -  x3uJn  +  Au;n2,  E12  =  el2  -  -x3uJi2  +  -A2tZ712 

1     _         1      _         1 

E\3    =     -x2w12  +  -A3u;i2  +  -\wnwi2 
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<f>2     =     ~^2wl2  +  -A3ti?i2  -  -Xwnwu 
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03     =     -~x3u;12  -  -A2it;12. 

Here  A(x2,  x3)  is  the  Saint- Venant  warping  function  for  the  beam  cross  section;  subscripts 
on  A  denote  partial  differentiation  with  respect  to  x2  and  x3.  Bars  over  a  symbol  denote  its 
value  at  the  beam  axis;  e.g.,  the  axial  displacement  of  the  beam  centroid  is  U(xi)  —  cuJi(xi). 
Substituting  these  relations  into  (1)  and  neglecting  higher  than  cubic  terms  in  the  dis- 
placements (these  are  not  needed  in  our  subsequent  analysis),  we  obtain  a  lengthy  expression 
for  the  potential  energy  which  forms  the  basis  for  our  subsequent  analysis.  Among  all  the 
functions  satisfying  the  geometric  or  natural  boundary  conditions  the  one  which  causes  the 
potential  energy  to  be  a  minimum  is  the  equilibrium  state.  We  suppose  that  the  outer  edges 
of  the  plate  are  free  to  displace  in  the  horizontal  plane  but  restrained  in  the  vertical  direction, 
so  the  geometric  boundary  conditions  are 

u;(0,x2)  =  u>(a,x2)  =  0  (2) 

Note  that  the  case  of  simply  supported  edges  is  obtained  by  setting  M  =  0. 


3.     PREBUCKLING  SOLUTION 

The  prebuckling  equilibrium  state  is  denoted  by  (u,io).    It  follows  from  assumptions 
(ii)-(v)  that  the  potential  energy  in  the  prebuckling  state  is 

,o   o,       ['    £(^  +  (6)2,      (6£>  +  £722)  £>„  o   c 

^"•"l-/,,    [ 2 + 2 ^B"*"  (3) 
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Here  722  is  the  moment  of  inertia  of  the  beam  section  about  the  x2-axis: 

-^22  =  /  /b«m  xl  dx2  dx3  . 
The  prebuckling  displacements  are  determined  by  the  variational  equation 

SP  =  0 

Taking  the  variation  of  (3),  and  integrating  by  parts  with  respect  to  £j,  we  obtain  the 
differential  equations 

(bD  +  EI22)  wuu  —EAc  ulu  —pb  =  0 
-EAc  wu  +E{A  +  tb)  ^  +a(A  +  tb)  =  0 

and  boundary  conditions 


(4) 


at  x  =  0  and  x  =  a  :    (bD  +  EI22)  Wn  -EAc  Ux  —a Ac  +  Mb  =  0  (5) 

The  solution  to  the  linear  boundary  value  problem  (2),  (4),  (5)  is 

o                 a       Acwu  ,c. 

Ui  =  -E  +  irnb  (6) 

o  px^(a  —  Xi)(a2 -\- axi  -  xl)  Mxx(a  —  xx) 

W    =     E A1? + E 7^~  (?) 

Note  that  (6)-(7)  reduce  in  the  case  I22  =  A  =  0  to  the  well-known  exact  solution  for  an 

isolated  wide  plate,  and  in  the  case  t  =  0  to  the  well-known  exact  solution  for  an  isolated 
beam. 


4.     BUCKLING  SOLUTION 

According  to  the  energy  criterion  of  elastic  stability,  the  prebuckling  equilibrium  state  is 
stable  if  and  only  if  the  energy  functional  which  represents  the  increase  of  the  total  potential 
energy  in  a  displacement  field  to  some  slightly  adjacent  state  (w  +w)  is  non-negative: 

P[w+w]-P[w]>0.  (8) 

Since  the  prebuckling  state  is  an  equilibrium  state,  the  terms  in  (8)  which  are  linear  in  the 
incremental  displacement  w  must  vanish.  It  follows  that  the  terms  Q[w]  in  (8)  which  are 
quadratic  in  the  incremental  displacement  must  be  non-negative: 

Q[w]  >  0 . 

The  critical  case  of  neutral  equilibrium  occurs  when  there  exists  a  buckling  mode  wcr  satis- 
fying 

QKr]  =  o  (9) 

Q[w^wa]>0.  (10) 

The  eigenvalues  crcr,pcr,  and  Mcr  which  render  (9)  zero  are  the  critical  buckling  loads. 

From  the  first  integral  in  (1)  and  assumption  (vi),  the  quadratic  functional  for  the  plate 

is 

atw\  "* 
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wt-,  win       ,  .     0 
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dx\dxi.         (11) 


2 
From  the  remaining  integrals  in  (1),  the  quadratic  functional  for  the  beam  is 

n       eh,  r_2       .{-EH3+ATTb  +  GH7)  r-  -^    AGJ-°i)  r-2  4 

Qb..m  =  — 2~y    ix;112dx1  + ^ J    wuw{idxl+ /    w^dx, 

(12) 
Here  /  is  the  polar  moment  of  inertia  about  the  zraxis,  J  is  the  Saint- Venant  torsion 
constant,  and  H\,  H2,  H3  are  constants  defined  by  the  following  integrals  over  the  beam 
cross  section: 


1    =  Lt.m(xl-\-xl)dx2dx3 

J      •'•fclion 

J      =       /Lm[(*2  +  A3)2  +  (*3-A2)2]Jx2<fx3 
J     •'tection 

Hi     =     /  Lm  A2</x2  <fx3 

H2    =     I  L>m  \xz{x\  +  x2  -  A2  -  A2)  +  2A(x2  +  A3)]<fx2  dx3 

J     •'•ection 

7/3    =     I  /  /»<-  ar3^2  "  Aa)2  +  (*3  +  A*)*  +  2(x2  +  x2  -  A2  -  A2)]</x2  <fx3 

The  total  quadratic  functional  for  the  plate  plus  the  beam  is  the  sum  of  (10)  and  (11). 

Next,  we  calculate  the  cross  section  properties  for  a  beam  composed  of  a  thin  web  and  a 
thin  bottom  and  top  flange.  The  web  has  thickness  tw  and  height  hw\  the  bottom  flange  has 
thickness  tb  and  width  hb\  the  top  flange  has  thickness  tj  and  width  hj.  The  Saint- Venant 
warping  function  for  this  thin-walled  cross  section  is: 

{x2{2hw  +  tj  -  x3)  flange 
x2X3  web 

— x2x3  bottom  flange 

Using  approximation  (viii)  we  obtain: 

A    =    tbh(,+  tu.hu,  +  t/hf 

Ac    =     -^— ^  +  tjhjhw 

twh3 
In    =     JyHi  +  if  h2w  hj 

.  hh\      twkl      tfh*f  2 

1    ~    T2"  +  "  +  -l2-~M/^ 

tjh       tlK       tjhj_ 
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tjhlh) 


Hx    = 


12 
t3jhwhj   ,  tlhl 


H2    =    ^-^  + 

Hz    m    h^  +  tfhlhf  +  ^ 


In  accordance  with  approximation  (vii)  and  the  boundary  conditions  (2)  and  (5),  we 
represent  the  incremental  buckling  displacement  by  the  following  shape  (an  arbitrary  multi- 
plicative constant  has  been  set  equal  to  1): 

w  =  sin sin  — —  ,    m=  1,2,3- ••  {io) 

a  b 

Substitution  of  this  buckling  mode  into  (11)  plus  (12)  and  application  of  the  inequality  (10) 
leads  finally  to: 

1orn      E„  A'<?\,        +  n      E,,  A2c\  +    2l  a    +  mb>   +        a>        +  °J 

„  <     12|g+i(/"-jTr6"      fl+i""-^ 

"  -  ,   a3 

(14) 

Here  m  is  taken  to  be  the  integer  which  gives  the  lowest  value  of  a^  in  (14).  Note  that 
(14)  reduces  in  the  case  t^  =  tw  =  tj  =  0  to  the  well-known  exact  solution  for  an  isolated 
plate,  in  the  case  pCT  =  M^  =  t  =  0  to  our  previous  exact  solution  for  an  isolated  beam  [see 
formula  (28)  of  Danielson  et  al  (1990)],  and  in  the  case  pCT  =  M„  =  0  to  formula  (87)  (with 
0pbe  =  (°e)™be  and  lzs2  +  r  as  Hx)  of  Adamchak  (1979). 

5.     CLAMPED  EDGES 

In  this  section  we  consider  the  clamped  case  when  the  bending  moment  M  is  not  pre- 
scribed, but  the  rotation  at  the  edges  is  completely  restrained.  Then  the  additional  geometric 
boundary  conditions  are 

wl{0,x2)  =  w1(a,x2)  =  0  (15) 

Note  that  in  this  case  each  material  particle  on  an  edge  cross-section  is  totally  restrained 
from  any  motion  at  the  buckling  point. 

The  prebuckling  solution  to  the  linear  boundary  value  problem  (2),  (4),  (15)  is  (6)  and 

*-    pt~xifAv  <i6>    i 

24'p  +  I^-T^)l 


Note  that  (6)  and  (16)  reduce  in  the  case  722  =  .4  =  0  to  the  well-known  exact  solution  for 
an  isolated  wide  plate,  and  in  the  case  t  =  0  to  the  well-known  exact  solution  for  an  isolated 
beam. 

In  accordance  with  approximation  (vii)  and  the  boundary  conditions  (2)  and  (15),  we 
now  represent  the  incremental  buckling  displacement  by  the  following  shape: 

w  =  (1  —  cos )  sin  — —  ,      T7i  =  2,4,0,- •  •  (i'j 

a  b 

Substitution  of  this  buckling  mode  into  (11)  plus  (12)  and  application  of  the  inequality  (10) 
leads  finally  to 


ac,< i ^+1^ jj, (18) 

Note  that  (18)  reduces  in  the  case  pCT  =  t  =  0  to  our  previous  exact  solution  for  an  isolated 
beam  [see  formula  (28)  of  Danielson  et  al  (1990)]. 


6.     IMPERFECTIONS 

In  this  section  we  suppose  that  the  structure  has  an  initial  normal  deflection  in  the  shape 
of  the  prebuckling  normal  displacement  caused  by  the  pressure  p.  Specifically,  when  M  is 
prescribed  on  the  edges  the  initial  normal  deflection  is 

16xi(a  —  X\)(a2  +  aii  —  x\) 


5a4 
while  for  clamped  edges  the  initial  normal  deflection  is 

16x5(0  -  xj)2 


W 


a4 
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We  also  assume  that  the  amplitude  W  of  the  initial  displacement  is  less  than  the  maxi- 
mum thickness  of  the  structure.  Then  the  prebuckling  displacements  are  still  given  by  our 


previous  solutions  (6)-(7)  or  (16),  and  the  quadratic  functional  Qpiate  is  still  given  by  (11). 

The  only  effect  of  this  prebuckling  deflection  is  to  create  a  new  term  in  Qbeam  which  is  the 

same  as  the  middle  integral  in  (12)  with  the  prebuckling  normal  displacement  replaced  by 

the  initial  displacement,  and  we  can  use  our  previous  calculations  to  evaluate  this  integral. 

We  thereby  find  that  the  effect  of  this  imperfection  is  to  add  an  additional  term  to  our 

previous  formula  for  a^: 

M  case: 

1  I  Ac 

32£(1  -  ~)(H3  -  i^j)W 

Ttlml A  -f  to 

aCT  =  aCT\W  =  0]  + ^73 (19) 
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Clamped  case: 


I  Ac 
c„  =  *„[W  =  0]  + *W (20) 

Note  that  for  an  asymmetrical  structure  torsional  deformation  of  the  stiffeners  may  ini- 
tiate upon  application  of  the  slightest  load,  so  bifurcation  may  not  be  able  to  be  used  as  the 
buckling  criteria  [see  Ostapenko  and  Yoo  (1988)]. 


7.     RESIDUAL  STRESS 

The  simplest  way  to  account  for  residual  stress  is  to  assume  that  the  plate  is  subject  to 
a  uniform  compressive  residual  stress  5,  while  the  beam  is  subject  to  a  counterbalancing 
distribution  of  residual  stress  ar(x2,x3)  [see  Hughes  (1983)].  The  only  effect  of  this  residual 
stress  is  to  create  a  new  term  in  QpUit  which  is  the  same  as  the  last  integral  in  (11)  with 
a  replaced  by  5.  Note  that  the  analogous  term  in  Qbeam  is  zero  because  w{x\)  =  0  for  the 
assumed  mode  shapes  (13)  and  (17): 


-III 


<jtw\ 

——axi  dx2  dx3  =  0 

beam         •- 
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We  thereby  find  that  the  effect  of  this  residual  stress  is  to  add  an  additional  term  to  our 
formulas  for  acr\ 

<rCT  =  <rcr[S  =  0] §-7-  (21) 

CONCLUSIONS 

Simple  analytical  formulas  that  include  the  effects  of  combined  loading,  various  boundary 

conditions,  imperfections,  and  residual  stress  do  not  appear  to  exist  in  the  literature.    As 

an  example  of  the  numerical  predictions  of  our  formulas,  let  us  assume  the  following  typical 

parameter  values  [taken  from  Smith  (1975)]: 

E  =  30,000  ksi  tb  =  0 

v  =  .3  hb  =  0 

a  =  48  in  tw  =  .28  in 

6  =  24  in  hu,  =  5.5  in 

t  =  .31  in  tj  =  .56  in 

hj  =  3.1  in 

For  this  example,  the  formulas  (14),  (19),  (21)  and  (18),  (20),  (21)  reduce  to 

M  case  (m  =  1): 

aCT  =  49  +  60/>cr  +  AbMcr  +  55W  -  .765  (22) 

Clamped  case  (m  =  2): 

acr  =  105  +  6.6Pcr  +  30VV  -  .765  (23) 

The  collapse  loads  of  a^  =  27.8  psi  (p^  =  0)  and  aCT  =  27.1  psi  (p^  =  .015  ksi)  measured 
by  Smith  (1975)  on  ship  grillages  may  be  accounted  for  by  choosing  appropriate  values  of 
McrtW,  and  5  in  (22)  or  (23). 

We  have  made  an  attempt  to  verify  the  accuracy  of  some  of  the  approximations  upon 
which  our  analysis  is  based.  For  instance,  in  the  M  case,  we  included  the  effect  of  beam 
cross-sectional  deformation  by  allowing  the  web  to  undergo  a  lateral  buckling  displacement 

v  =  sin (^1X3  -f  C2X3)  , 

a 
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where  C\  and  Ci  are  constants  determined  by  minimizing  the  critical  axial  stress  a^.  For 
the  parameter  values  listed  above,  this  effect  turned  out  to  be  negligible.  It  is  possible  to 
invent  unusual  cases  in  which  the  cross-sectional  deformation  is  of  importance,  but  for  most 
practical  dimensions  the  assumption  of  a  rigid  cross-section  seems  ok. 

For  another  instance,  in  the  clamped  case,  we  added  together  the  m  =  2  and  m  =  4 
normal  buckling  displacements  (17): 

.  27TX!  4ttxi       .    nx2 

w  =   1  —  cos (-  63(1  —  cos )J  sin  — —  , 

a  a  b 

where  C3  is  a  constant  determined  by  minimizing  the  critical  axial  stress  aCT.  For  the 
parameter  values  listed  above,  this  effect  also  turned  out  to  be  negligible.  For  most  practical 
dimensions  the  assumption  of  a  simple  buckling  mode  seems  ok. 

At  any  rate,  if  we  include  the  additional  displacement  functions  needed  for  the  above 
effects,  it  doesn't  seem  possible  to  obtain  simple  formulas  for  the  unknown  coefficients  or 
buckling  loads. 
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